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ABSTRACT 

Aims. Recurrent nova systems like RS Oph have been proposed as a possible channel to Type la Supernova explosions, based on the 
high mass of the accreting white dwarf. Additional support to this hypothesis has been recently provided by the detection of circum- 
stellar material around SN 2006X and SN 20071e, showing a structure compatible with that expected for recurrent nova outbursts. 
In this paper we investigate the circumstellar environment of RS Oph and its structure, with the aim of establishing a firmer and 
independent link between this class of objects and Type la SN progenitors. 

Methods. We study the time evolution of Can, Nai and Ki absorption features in RS Oph, before, during, and after the last outburst, 
using multi-epoch, high-resolution spectroscopy, and applying the same method adopted for SN 2006X and SN 20071e. 
Results. A number of components, blue-shifted with respect to the systemic velocity of RS Oph, are detected. In particular, one fea- 
ture strongly weakens in the first two weeks after the outburst, simultaneously with the disappearance of very narrow P-Cyg profiles 
overimposed on the much wider nova emission lines of H, He, Fe ii and other elements. 

Conclusions. We interpret these facts as the signature of density enhancements in the circumstellar material, suggesting that the 
recurrent eruptions might indeed create complex structures within the material lost by the donor star This establishes a strong link 
between RS Oph and the progenitor system of the Type la SN 2006X, for which similar features have been detected. 
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1. Introduction 



RS Ophiuchi is a symbiotic recurrent nova (Kenyon'1986'), that 
has undergone a number of recorded eruptions, with an inter- 
outburst period of 10-20 years. The binary system is composed 
by a K4-M0 red giant (RG) star in a 456 day orbit with a 
white dwarf (WD) of mass close to the Chandrasekhar limit 
(Dobrzycka & Kenyon 1994; Shore et al. [T9961 Fekel et al. 
12000; Brandi et al. ,2009) . For this reason, RS Oph has been 
proposed as the prototype for one of the possible channels to 
Type la Supernova explosions (see for example Hachisu & Kato 
2000 2001 ; Hernanz & Jose 2008 ). The WD is supposed to ac- 
crete part of the material lost by the RG via stellar wind, until it 
reaches the conditions for thermonuclear burning at its surface. 
When this happens, the WD ejects 10 to 10"* M© at veloci- 
ties of -5,000 km s-i (Hachisu & Kato 2001 ; Yaron et al. l2005T l. 
producing a shock wave that propagates into the pre-existing RG 
wind and gives rise to strong X-ray emission, as observed in the 
1985 eruption (Mason et al. 1987J . The last outburst of RS Oph 
took place on 2006 February 12.83 UT (Narumi et al. l2006l) . and 
prompt radio observations have shown that the nova ejection was 
bipolar (O'Brien et al. 2006). Besides confirming that the blast 
wave clearly deviated from spherical symmetry, the X-Ray data 
(Bode et al. 120061: Sokoloski et al. I2006I I have indicated a low 

Send offprint requests to: F. Patat 
* Based on observations obtained at ESO-La Silla. 



mass for the ejected material which, in turn, favors a high mass 
for the WD (~1.4 Mq), giving further support to the proposed 
connection between RS Oph and Type la progenitors (Sokoloski 
et al. 2006 ). One important question that remains unanswered, is 
whether the hot component of RS Oph is a O-Ne-Mg WD rather 
than a C-O WD, for there are no direct observations that could 
help to distinguish between the two. In this case, the WD would 
not experience a thermonuclear explosion, leading to a SN la, 
but rather undergo an electron-capture core collapse, producing 
a neutron-star remnant (Nomoto & Kondo [T99Tj . when its mass 
gets close to the Chandrasekhar mass. 



In addition to the clear detection of circumstellar matter 
(CSM) very close to the system (O'Brien et al. 2006; Sokoloski 
et al. |2006] |. lijima (l2008l |2009J has reported the presence of 
a narrow Na i D absorption component, which is blue-shifted by 
about 33 km s ' with respect to the systemic velocity of RS Oph. 
He proposed that this feature might arise in the wind of the RG, 
hence arguing for its circumstellar nature. Stimulated by this 
finding and by the detection of Na i D time-variant features in the 
Type la SN 2006X (Patat et al.l2007] l. SN 20071e (Simon et al. 
20091 ) 1999 cl (Blondin et al.'2008^, and SN 2006dd (Stritzinger 
et al. 12010! ) we performed a thorough analysis using pre/post- 
outburst multi-epoch, high-resolution spectroscopy of RS Oph, 
whose results we are reporting here. 
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Table 1, Log of RS Oph FEROS observations. 



Date MJD 


Phase 




Epoch 


Airm. Exp. Bar.Corr. 






(km s-') 


(days) 




(s) 


(km s-') 


2004-04-11 53106.39 


0.48 


-16.6 


-672.4 


1.08 


1200 


+26.1 


2004-06-05 53162.23 


0.60 


-13.5 


-616.6 


1.09 


1200 


+6.3 


2004-06-06 53163.29 


0.60 


-13.3 


-615.5 


1.15 


1200 


+5.7 


2004-08-31 53248.05 


0.79 


+4.1 


-530.8 


1.13 


1200 


-26.9 


2006-02-14 53780.37 


0.96 


-1-16.1 


+ 1.54 


2.09 


120 


+24.7 


2006-02-25 53791.37 


0.98 


+ 16.6 


+ 12.54 


1.57 


30 


+27.1 


2006-03-07 53801.40 


0.00 


+ 16.7 


+22.57 


1.22 


90 


+28.4 


2006-03-18 53812.42 


0.03 


+ 16.5 


+35.59 


1.10 


180 


+28.7 


2006-03-31 53825.41 


0.06 


+ 15.7 


+46.58 


1.08 


300 


+27.9 


2008-02-24 54520.35 


0.58 


-14.5 


+741.52 


1.72 


4840 


+26.9 


2008-02-25 54521.35 


0.58 


-14.5 


+742.52 


1.75 


4200 


+27.0 


2008-02-26 54522.36 


0.58 


-14.5 


+743.53 


1.78 


4650 


+27.2 



Note: Epochs are computed from the 2006 outburst (February 12.83 UT, 
Narumi et al. i2006i) . Phase and radial velocities for the M giant (r^) are 
computed using the ephemeris by Fekel et al. C2000) . 



2. Observations and Data Reduction 

RS Oph was observed in the pre-outburst, outburst and post- 
outburst phases using the Fiber-fed Extended Range Optical 
Spectrograph (FEROS; Kaufer et al. [T999T l mounted at the ESO- 
MPG 2.2m telescope. The pre-outburst data set was presented 
by Zamanov et al. (2005 ), while the outburst data are going to be 
thoroughly discussed in a forthcoming paper, to which we refer 
the reader for a full report on the high-resolution spectroscopic 
follow-up of RS Oph (Mason et al. in preparation). In order to 
cover the quiescence phase following the outburst, an additional 
set of data was obtained in February 2008, about two years after 
the 2006 eruption. 

The instrument delivers a resolving power A/AA ~ 48,000, 
covering the spectral range 3700-9200 A. The data were pro- 
cessed by the FEROS pipeline, which includes the correction for 
Earth's motion. The first set of data was obtained on four epochs 
between 672 and 530 days before the February 2006 outburst, 
whilst the second covers 5 epochs between 1.5 and 46.6 days af- 
ter the eruption. Finally, RS Oph was observed on three consec- 
utive nights about 740 days after the eruption, when the system 
was well back into the quiescence phase, and the mass accretion 
had resumed (Worters et al. l2007l . The spectra obtained on these 
three nights were stacked in order to increase the signal-to-noise 
ratio. The resulting spectrum is presented and discussed here for 
the first time. 

A journal of the observations is given in Table [T] The pres- 
ence and position of telluric lines, produced by H2O and O2 ab- 
sorption bands, was checked using a synthetic spectrum com- 
puted by the Reference Forward ModelM Gas column densities 
were estimated using VPGUESS0 and refined with WPFYf\ 



3. Inter- and/or circumstellar lines towards RS Oph 

Pronounced Ca 11 H&K, Na i D, K i /17665, and K i ^7699 are de- 
tected in the high-resolution spectra of RS Oph, as it is shown in 
Fig. [T] While the K i A1699 profiles are rather severely affected 
by telluric absorptions, these are practically negligible for the 
other lines. Besides two intense features at heliocentric velocity 



' http://www.atm.ox.ac.uk/RFM/ 
2 VPGUESS was developed by J. 



-o 

0) 




Liske and can be downloaded at 
jhttp://www.eso.org/~jliske/vpguess/index.html 

^ VPFIT was developed by R.F. Carswell and can be downloaded at 
|http://www.ast.cam.ac.uk/~rfc/vpfit.html| 



-100 -50 

Observed Heliocentric Velocity (km s"^) 

Fig. 1. Interstellar lines along the line-of-sight to RS Oph: Can 
K (upper), Na I Di (middle), and K i /17665 (lower) on day -672. 
The vertical dashed lines mark the main components (see text). 
The dashed-dotted line indicates the systemic velocity (-40.2 
km s Fekel et al. l2000l l. The green curve is a synthetic telluric 
absorption spectrum matched to the observed spectrum. 



-12 km s ' (#4) and +2 km s ' (#5), arising in the local spiral 
arm (Hjellming et al. 1 19861 1, the data reveal at least three distinct 
component^at -77 (#1), -63 (#2) and -46 km s"' (#3), respec- 
tively. The same components, though with different intensities, 
are visible in all Can, Nai, and Ki features (Fig. [TJ. Besides 
those already mentioned, other shallower features are most likely 
present between -40 and -20 km s"', but they contribute only 
marginally to the overall profile. 

As the systemic velocity of RS Oph is -40.2 km s"' (Fekel 
et al. 120001 see also Brandi et al. I2009I I. components #1 to #3 
must arise in material coming towards the observer, at velocities 
where no interstellar component is expected (KeiT & Westerhout 
1965i s ee also Cassatella et al. ir9B5i Hjellming et al. |1986t 
Snijders[l987]l. 

This is a most intriguing fact, especially if it can be proven 
that these absorptions arise in the circumstellar environment. In 
this context, the first question to be answered is about the possi- 
ble/probable presence of a photospheric absorption component, 
which may further complicate the picture. 



4. The spectrum of the giant during quiescence 

The spectral type of the donor star in the RS Oph system, based 
on the study of molecular bands in the near-IR, is uncertain and 
ranges from K4 III to MO III (Murset & Schmid 1999) . The ef- 
fective temperature is estimated to be in the range 4100-4400 K, 
and the metallicity is [M/H] =0.17+0.1. Additionally, the star is 
reported to show an anomalously high Li abundance (Wallerstein 



* In the spectra presented by lijima ( 120081 12009b . obtained on 
February 18.2 (+5.4) and 19.2 (+6.4) with a resolution A/AA -10,000, 
these components are severely blended. 
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Fig. 2. Time evolution for Na I Di (lower panel) and Cai /16717.7 
(upper panel) on three quiescence epochs (-672, -530 and 
+742). The vertical grey strip indicates the radial velocity range 
of the M giant during the orbital period (±16.2 km s Fekel et 
al. |2000). The dashed curve in the lower panel is the line pro- 
file of the K giant star HD32887 (Da Silva et al. |2006] |. veiled 
to 30% of its original depth and placed at the velocity deduced 
from Cai i6718 for day -530. 

et al. 120081 see also the discussion in Hernanz & Jose 120081 and 
Brandietal.5009). 

Because of the presence of a hot component, the spectrum 
of the giant is heavily veiled (Bruch 1986), so that below 4000 
A this is dominated by the hot component, while at 5500 A 
the contribution of the cool component is between 30 and 50% 
(Dobrzycka et al. 1996). For these reasons it is natural to expect 
the strong and wide photospheric Na i D lines typical of a cool 
giant to be visible in the composite spectrurrj^ 

Therefore, a possible explanation for the changes in the Na i 
D absorption features during the quiescent phase (see Fig. |2] 
lower panel) is the variation of the lines in the M giant spec- 
trum. These might be explained in terms of wavelength shifts 
due to the orbital motion, but also as intrinsic line variability. 
In fact, the hot components of symbiotic systems are known to 
heat the facing hemispheres of the red companion, hence pro- 
ducing temperature changes in their atmospheres large enough 
to change the lines shapes (see Kenyon et al. ll993l and references 
therein). Indeed, in the case of RS Oph both mechanisms might 
be at work. This is shown in the upper panel of Fig. |2] where 
we plotted the region of the Cai /16718 line for the first (-672 
days) and the last (-530 days) pre-outburst epochs in our data 
set. Besides the velocity shift observed between the two epochs 
(~25 km s"'), the Ca I line shows a significant depth increase, 
a behavior that is common also to other stellar lines (like Li i 

^ Marked Na i D lines, supposedly arising from the cool component, 
were reported by Bruch ('1986'), with a global equivalent width of about 
3 A. Nevertheless, given the low resolution of Bruch's data, most of 
the absorption he measured is most probably caused by the interstellar 
component. 



Fig. 3. Comparison between the evolution of Can 8662.14 (red) 
and Nai D] (black) during three quiescence epochs. The blue 
curve traces the Cai /16718 line, while the blue arrow marks the 
estimated M giant radial velocity. 

/16707, very neatly detected in our spectra). This might explain, 
at least partially, the evolution of the Na i D absorption profiles 
observed during the quiescence phase (Fig.|2] lower panel). 

An important fact to be noticed is that Na i D and Ca ii H&K 
regions are most probably aff'ected by the presence of compo- 
nents arising in the circumbinary environment, which are char- 
acterized by the complex profiles seen in other prominent lines 
like H and He (see for instance Zamanov et al. 12005 1 Brandi 
et al. 12009 l l. Because of the presence of interstellar and pos- 
sible circumstellar components, it is very difficult to distin- 
guish these components in Nai D, Can H&K, and Ki lines. 
Nevertheless, this becomes feasible for other strong lines, where 
the circum/interstellar components are absent, like the near-IR 
Can triplet. An example is presented in Fig. [3] where the evo- 
lution of Can /18662.14 is compared to the one of Nai Di. The 
changes seen in the NIR-Ca ii lines are followed by those in the 
Na I profile and can be explained, at least partially, by the radial 
velocity shift of the M giant, traced by the Cai /16718 line. 

Therefore, during quiescence there certainly is a stellar com- 
ponent, possibly characterized by a complex structure, which af- 
fects the observed profile of the narrow absorption lines arising 
in the inter/circumstellar environment. 

5. Evolution of absorption features during outburst 

The time evolution during one pre-outburst (-672 days), two 
outburst (+1.5, +12.5 days), and one post-outburst (+742 days) 
epochs is presented in Figs. [4] [5] |6] for Can H, Nai Di, and Ki, 
respectively. The most striking fact shown by these plots is the 
marked evolution of component #3. 

For the specific case of Na i, in the pre-outburst epoch this 
component is almost completely saturated, while on day +1.5 
its depth has decreased by about a factor 2. The line appears 
to have weakened further on day +12.5, after which it does not 
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Fig. 4. Evolution of the Can H profile on days -672, +1.5, 
+ 12.5 and +740. The horizontal bar indicates the radial veloc- 
ity range of the M giant during a full orbit. 



show any significant evolution up to the last epoch covered by 
our data (+46.6). The column density retrieved from component 
#3 changes from logA^ ~12.9 on day -672 to ~12.1 on day 
+ 1.5, to finally reach ~11.3 on day +22.6, implying a variation 
of over a factor 30. The line width (/7=FWHM/1.665) derived 
for this component ranges between 3 and 4 km s"' during the 
post-outburst epochs, while for the pre-outburst epoch the best 
fit gives b ^8.5 km s This is a further indication that an addi- 
tional component is present prior to outburst, most likely due to 
the M giant intrinsic absorption (see Sect.|4]i. 

In general, the behavior seen in Can H&K is very similar 
to that of Nai D, with component #3 clearly decreasing in in- 
tensity, both between the pre-outburst (-672) and the first out- 
burst (+1.5) epochs, and between the first two outburst epochs 
(+1.5, +12.5). Importantly, the pronounced evolution displayed 
during the first two weeks after the outburst clearly demonstrates 
that this cannot be attributed solely to the disappearance of the 
photospheric component. In fact, during these phases the over- 
all luminosity of the system is ~100 times larger than that of 
the M giant, making its relative contribution to the total spec- 
trum negligible. This is illustrated in Fig. [7] where we plotted 
the comparison between the evolution of the near-IR Ca ii /18662 
and Nai Di across several of the epochs covered by our obser- 
vations. While during quiescence the complex profile of Can is 
clearly detected, this is completely missing during the outburst. 
On the contrary, barring the variation seen in component #3, the 
blue-shifted Nai components are always visible, ot all epochs. 
The same is true for the Can H&K features. 

Therefore, absorption components #1, #2 and #3 certainly 
arise outside of the stellar atmosphere, and outside the nova 
ejecta. In addition, the variation seen in component #3 tells us 
that at least that feature cannot be interstellar In this respect, 
an even more conclusive argument about the nature of all blue- 
shifted features comes from a closer inspection of Fig.|6] At vari- 
ance with what is seen for Na i D and Ca n H&K, the K i absorp- 
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Fig. 5. Same as Fig.^for Nai Di. 



tions disappear almost completely during the outburst (compo- 
nent #2 is still visible on days +1.5 and +12.5. See also Sect.|6]l. 
This definitely rules out an interstellar origin for these lines. 

A hint to the physical mechanism responsible for the ob- 
served line variations comes from the following consideration. 
The variations appear to be stronger for species with lower ion- 
ization potentials: the largest changes are seen for Ki (4.3 eV), 
while the weakest ones are observed for Can (11.9 eV), with 
Nai (5.1 eV) showing an intermediate behavior. This suggests 
that the weakening of absorption features is related to a change 
in the ionization balance, induced by the radiation field produced 
by the nova and/or by the interaction between the nova ejecta and 
pre-existing, circumstellar material. 

At variance with the SN case, where the system is completely 
disrupted by the explosion, the recurrent nova case offers the 
possibility of studying the CS environment also after the out- 
burst. With the aid of the +742 epoch, obtained when the system 
is supposed to be back to the quiescence phase (Worters et al. 
120071) . one can directly verify whether the outburst has modified 
the CSM in a sensible way. In this respect, the evolution shown 
in Figs.|4j|5] and|6]is revealing. Although the post-outburst epoch 
was obtained at a similar orbital phase of the pre-outburst epoch 
-672 (see Table [TJ, the profile is markedly different, with an en- 
hanced absorption in the blue. This is common to Ca n, Na i, and 
K I, pointing to a global modification of the circumstellar envi- 
ronment produced by the nova ejecta. 

This conclusion is strengthened by the analysis of the Ha 
profile presented in Fig.|9] The post-outburst profile departs from 
that seen during the 3 pre-outburst epochs, in that the absorption 
extends far more into the blue, reaching a systemic velocity of 
about -60 km s ' . Given the evolution shown in the pre-outburst 
phases, it is evident that the variation seen on day +742 cannot 
be explained in terms of the pure orbital motion of the M giant. 
The question as to whether this is a transient phenomenon or a 
more stable perturbation of the circumstellar configuration will 
have to await for follow-up high resolution spectroscopy. 
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Fig. 6. Same as Fig. [4] for Ki /17665. The absorption at — 1-65 
km s"' is a telluric feature. 



A very remarkable fact observed in conjunction with the 
weakening of component #3 is the disappearance of a narrow 
P-Cygni profile which is well visible in H, He, Oi and basi- 
cally all other broad emission lines on day +1.5. For example, 
on this epoch. Ha shows a rather complex, multi-component 
profile, with a FWHM of about 3400 km s"', indicating excep- 
tionally high expansion velocity for the nova ejecta. On top of 
this emission, a narrow P-Cyg profile is present, with the mini- 
mum of the absorption trough placed at a velocity of about -10 
km S-' (Fig. [s] upper panel). This is notably similar to that of 
the time variable Caii/Nai/Ki feature #3 (see also Fig. [8] lower 
panel). On the next epoch (day +12.5), this narrow feature has 
completely disappeared, leaving the pure underlying Ha broad 
emission profile with FWHM~1250 km s ' . 

The disappearance of the narrow P-Cyg absorptions soon af- 
ter the outburst indicates that the gas in which they originate 
must be close to the eruption site, in order to be influenced 
by the radiation field and/or by the direct interaction with the 
nova ejecta, that is known to start immediately after the outburst 
(O'Brien et al. ,2006. Sokoloski et al. .2006< . 

The most natural conclusion one can draw from all these 
facts is that the narrow, blue-shifted absorptions observed in RS 
Oph arise in circumstellar material. 



6. Discussion 

The time variations of the Na i, Ca 11 and K i lines observed in RS 
Oph are rather dramatic, and cannot be compared to the known, 
long timescale fluctuations observed along the line of sight to 
galactic stars (Crawford 1 1992l l. They are much more resemblant 
to those observed in the Type la SN 2006X (Patat et al. 2007 ) and 
20071e (Simon et al. 12009). These were interpreted as arising 
in the circumstellar material, whose physical conditions where 
probably modified by the explosion, both by the radiation field 
and by the direct interaction with the SN ejecta (see the discus- 
sion in Chugai 120081 for an alternative scenario. See also Patat 



et al. 120101) . To the best of our knowledge, short timescale vari- 
ations of narrow absorptions are known to this level of detail 
only for another object, i.e. the core-collapse, type Iln SN 1998S 
(Bowen et al. 2000 ). High resolution spectroscopy of this object 
has revealed a number of Nai D components within the host 
galaxy, the bluest of which deepened significantly during the 19 
days spanned by the two epochs available, implying a Na i col- 
umn density increase of >1 dex. SN 1998S has shown definite 
signs of ejecta-wind interaction in the form of narrow H and 
He P-Cyg profile lines and, remarkably, these features were de- 
tected at the same velocity of the Na i D variable component, and 
showed the same deepening with time. Bowen et al. (2000) con- 
cluded that the narrow and variable absorption is a signature of 
the outflows from the super-giant progenitor of SN 1998S, arising 
in a dense shell, expanding at about 50 km s"' (see also Fassia 
et al. 200 ll 

Given all the evidences collected in this work (see Sect. [5]), 
the most plausible possibility is that the time variant Na i, Ca 11, 
and K i features are produced within the circumstellar environ- 
ment of RS Oph. Since clear signs of interaction between the 
CSM and the nova ejecta were detected immediately after the 
eruption (O'Brien et al. 120061 Sokoloski et al. 2006), it is rea- 
sonable to argue that the weakening of component #3 is due to 
the interaction with the fast moving nova ejecta, which sweep 
away the pre-existing material. Under this hypothesis and given 
the fact that the interaction must have akeady started by the time 
of our second outburst epoch, an upper limit to the distance of 
this gas can be obtained from the ejecta velocity (~4,000 km s"') 
and the time of interaction (<12.5 days), and it turns out to be 
<4xl0''* cm. 

The velocity measured for the time-variable component #3 
(6.3 km s"') is significantly smaller than that expected for a RG 
wind. Even taking into account the combination of wind veloc- 
ity, rotational velocity, orbital velocity and hydrodynamical ef- 
fects, the radial speed of the outflowing material is expected to 
be between 20 and 50 km s depending on the viewing an- 
gle (Walder, Folini & Shore I2008T ). Therefore, it is plausible that 
component #3 arises within the circumbinary environment, at 
variance with components #1 and #2, which most likely arise 
much farther out (see the discussion below). 

The problem with this scenario is the absence of a P-Cyg 
profile in the Na, Ca and K lines, contrary to what is seen in H, 
He I, O I, Fe 11 and Si 11 lines. The most viable explanation is that 
Na I, Ca II, and K i lines arise in an asymmetric density enhance- 
ment lying along the line of sight, while H and other lines are 
generated both in an asymmetric and in a more symmetric mate- 
rial, possibly identified as a circumbinary equatorial wind. The 
reason why we do not see Na i, Ca 11, and K i lines arising from 
the symmetric component is plausibly a strong ionization due to 
their low ionization potential. In fact, all the ions for which P- 
Cyg profiles are observed have higher ionization potentials (H: 
13.6eV;Hei:24.6eV;Oi: 13.6 eV; Si 11: 16.3eV;Feii: 16.2 eV) 
than Nai (5.1 eV), Can (11.9 eV) and Ki (4.3 eV) and higher 
(by at least 1 dex) abundance than Na, Ca and K. Therefore the 
latter three species, being strongly ionized in the wind, are prac- 
tically unobservable. Only in the dense asymmetric component 
the optical depth in these absorptions is sufficient to produce de- 
tectable features. In this respect we note that the emergence of 
asymmetric, dense structures lacking point-symmetry, is a nat- 
ural consequence of the orbital motion of a mass-loosing red 
giant (Mastrodemos & Morris I1999D . In general, when dealing 
with a system like RS Oph, one should bear in mind that the 
dynamic conditions of the circumbinary gas are most likely non- 
stationary, because the flow is strongly affected by the orbital 
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Fig. 7. Comparison between the evolution of Can 8662.14 A 
(red) and Nai Dj (black) on days -672, -615, -531, +1.5, 
+ 12.5, and +742. The dashed vertical lines mark the expected 
radial velocity of the M giant (Fekel et al. i2000i) . 



Fig. 8. Comparison between the evolution of Ha (upper panel) 
and Na I Di (lower panel) on the first two, post-outburst epochs. 
The small insert in the upper panel shows the overall profile of 
the Ha line. The vertical dashed line marks the position of the 
time variable Nai D component #3. 



motion. This may result in the formation of transient absorbing 
structures in the line of sight at similar orbital phases, which 
make very difficult a proper modeling of the geometry and gas 
flow of the CSM. Hydrodynamical 3D-simulations by Walder et 
al. (2008) show that during quiescence the circumbinary matter 
of RS Oph can be described as an oblate spheroid on the scale 
of ~ 10'^ cm, while on smaller scales (~ 10'^ cm) the geometry 
is rather complicated, and essentially lacks point-symmetry. 

As we mentioned above, the early disappearance of compo- 
nent #3 is strongly suggestive of a circumbinary nature. In this 
respect, it is interesting to note that the so called cF-absorption 
system was destroyed during the outburst (Brandi et al. |2009l l. It 
then reappeared (accompanied by flickering) about 8 months af- 
ter, indicating the resumption of accretion (Worters et al. 2007). 
The cF-absorptions are believed to be related to the hot com- 
ponent of the binary system (Mikolajewska & Kenyon I1992I I. 
In addition, their measured velocities are lower than 10 km s"^ 
(Brandi et al. 2009), which is consistent with the low velocity 
deduced for component #3 (~6 km s"'). This establishes a fur- 
ther and independent link between this feature and the circumbi- 
nary material. Brandi et al. associate the cF-absorption system to 
the material streaming towards the hot component. We empha- 
size that if component #3 arises within this gas, that conclusion 
does not hold, for this feature is still visible on day +1.5, and it 
must therefore lie farther out. Indeed, at this age the radius of the 
nova shell is a!4.5xl0'^cm, i.e. larger than the binary separation 
(s;2xl0i^^ cm). 

From the spatially resolved HST observations of the 2006 
eruption, it is clear that the nova remnant has a bipolar morphol- 
ogy, resulting from the interaction of the outburst ejecta with 
CSM that is significantly denser in the equatorial regions of the 
binary than at the poles (Bode et al. 120071 ). In general, the inter- 
action of the nova ejecta with the CSM should essentially modify 
the density and velocity field of the CS environment formed dur- 



ing quiescence. The accelerated CS gas presumably acquires a 
broad range of densities and velocities, depending on the polar 
angle. In this respect, it is tempting to explain component #1 as 
the outcome of the interaction between the nova ejecta and the 
RG wind lost during the last quiescent phase. 

To provide a first order insight into this possibility, we con- 
sider a spherical approximation, and assume that the mass M^. 
is lost by the RG with a wind velocity m,j, in the course of the 
last inter-outburst phase (1986-2006). The nova ejecta, charac- 
terized by mass and kinetic energy E,,, sweep up the material 
lost during quiescence, and forms a shell with mass Mg + M„. 
This expands at a final velocity u, which is dictated by momen- 
tum conservation. To estimate this velocity one needs to express 
the initial momentum carried by the ejecta as a function of 
and Mg. To this end, we assume that the nova ejecta start to 
expand homologously, with an exponential density distribution 
p oc exp (-v/v()), where the scale velocity vq is determined by Eg 
and Mg. Adopting typical values Mg - 10"^ Mq and Eg - 10"*^ 
erg, one gets vo = 9.1 x 10'' cm s"' (note that the outer ejecta of 
course have much larger velocities). The ejecta momentum, inte- 
grated over the sphere for the adopted ejecta structure, is IMgVQ. 
Momentum conservation then leads to the following relation be- 
tween mass ratio fi = Mg/M„, and velocities 

u - u„ 
2V() - u 

Adopting a typical RG wind velocity m„, = 20 km s \ a final 
velocity of the swept-up shell m = 50 km s"', the ejecta parame- 
ters Eg = 10'^^ erg and Mg = 10"^ Mq, from the above expression 
one obtains Mg/M„ = 0.017, or M,y » 6 x lO"*" Mq. The latter 
implies a RG mass-loss rate Mrg=M,v/20 yr » 3 x 10"^ Mq yr"' , 
which is quite a sensible value for a RG in a symbiotic system 
(Seaquist & Tavlor [T990t . 
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Fig.9. Evolution of the Ha profile on days -672, -615, -531, 
+ 1.5, and +742 (from top to bottom). The dashed vertical lines 
mark the expected radial velocity of the M giant (Fekel et al. 
I2OOO1) . 

Despite our crude approximations, this estimate indicates 
that component #1 might be related to the wind material acceler- 
ated by the nova ejecta. If the absorption is produced by the shell 
formed by the previous 1986 outburst, then the absorber with a 
velocity of ~50 km s"' is located at a radius ~3xlO'^ cm. In this 
scenario, this spatial scale would also be equal to the distance Ar 
between subsequent shells. 

This simplified picture for component #1 neglects two im- 
portant issues: i) the deviation of the CSM from spherical sym- 
metry, which results in the dependence of final velocities on the 
polar angle; ii) the Rayleigh-Taylor instability of the swept-up 
shell, which results in shell fragmentation. The latter is essential 
for understanding the absence of scattered emission in compo- 
nent #1 (see above). Indeed, if the covering factor of fragments is 
C ^ 1, then the intensity of the scattered emission is suppressed 
by the same factor C with respect to that arising in the spherical, 
smooth shell. On the other hand, the absorption can be rather 
strong (i.e. comparable with that of the spherical component) if 
the absorbing cloud accidentally gets into the line of sight. 

In fact, the absorption can be produced also by more dis- 
tant CS clouds formed by previous outbursts and mass-loss. The 
Ca II column density required to produce a Ca n /13968 line with 
an optical depth t-1 is 2.8 x 10'^ cm"^ (for b - 10 km s"'). If 
calcium is mostly Ca 11, then for a solar abundance this line can 
be detected, provided the cloud produced by the shell fragmenta- 
tion lies at a distance of about 2xlO"'C"'^^ cm. For instance, in 
the case of C ~ 10"^, Can would be detectable up to a maximum 
radius r,„ax ~ 2x lO'^ cm. 

Based on the available data, lijima (120081 12009I I identified 
the blue-shifted component at systemic velocity -33 km s"' as 
possibly arising within the wind of the RCj^ Our data, which 

Note that lijima l|2008l|2009ll used a systemic velocit y of -35.8 km 
s ' for RS Oph. Adopting the value of Fekel et al. ( 12000b . the systemic 
velocity of the blended Nai component becomes -29 km s~'. 




Relative Velocity (km s"') 

Fig. 10. Comparison between high-resolution spectroscopy of 
three Type la SNe and RS Oph in the region of Nai D2. For 
presentation the velocity scales refer to the center of the deep ab- 
sorption, most likely associated to the local spiral arm. The data 
of SN 1991T, 1998es and 2006X are from Patat et al. (l2007] l.The 
light colored curve in the bottom panel is an atmospheric model. 



have a spectral resolution almost a factor 5 higher, show that 
multiple components are present. As we have been arguing 
throughout this paper, this is strongly suggestive of a more com- 
plex configuration, in which the interaction between the recur- 
rent outbursts and the RG wind plays a fundamental role in shap- 
ing the CSM. Although we cannot make a firm statement, it is 
possible that component #2, characterized by a systemic velocity 
of about -23 km s ', is the signature of the RG wind. Another 
possibility is that it is the left-over of a previous outburst, sim- 
ilarly to what we proposed for component #1. At variance with 
components #1 and #3, #2 is still present in the Ki profile after 
the outburst (see Fig. [6]). This indicates that the gas responsible 
for this absorption is placed at a sufficiently large distance that 
its ionization balance is not too much affected by the eruption. 

The presence of component #3 after the nova outburst can 
be explained by the resumption of the accretion, so that its re- 
appearance in our latest epoch is not inconsistent with it being 
formed within the circumbinary material (r <10''* cm). On the 
other hand, the presence of #1 and #2 at the same epoch is hard 
to explain if they arise too close to the WD. Since the nova shell 
reaches a radius of a;10'^ cm around one month after the erup- 
tion, they would have been swept away rather soon, had they 
been placed at shorter distances. This suggests that the absorb- 
ing material is placed at a sufficiently large radius that the ejecta 
are decelerated well before reaching it. In the light of these con- 
siderations, and given the persistence of the Ca 11 H&K and Na i 
D lines during the outburst, the most probable cause for the vari- 
ability of components #1 and #2 is photoionization. 

If this is correct, then the recurrent outbursts do not destroy 
the slow-moving shells produced by the previous events, imply- 
ing that systems like RS Oph are able to produce long-lasting 
structures in their circumstellar environment. Structures that sur- 
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vive long enough to show up during the most definitive and 
extreme event, when the WD is finally incinerated by the ther- 
monuclear runaway. 

7. Conclusions 

Since the introduction of the mass-accreting WD scenario to ex- 
plain Type la SN explosions (Whelan & Iben .l973j . the exact 
nature of the prog enitor system/s has been extensively debated 
(see Branch et al. ^T995'i Livio '2000', Parthasarathy et al. l2007l 
Tutukov & Fedorova 2007 for extensive reviews). Among the 
various possibilities, symbiotic systems have been proposed as a 
viable channel (Munari & Renzini il992i ). In particular, because 
of the mass of the WD hosted by RS Oph, which is close to the 
Chandrasekhar limit, this system has been proposed as a strong 
candidate (Hachisu & Kato 2000 2001 ). The detection of CSM 
material around SN 2006X and 20071e has strengthened this 
relation, showing the presence of multiple, high-density, low- 
velocity clumps of neutral material along the line of sight (Patat 
et al. 122)7; Simon et al. 120091 1. 

The analysis presented in this paper, stimulated by the find- 
ings reported by lijima ( 12008! |2009), leads us along the same 
lines. If the presence of narrow, blue-shifted absorptions at ve- 
locities where no interstellar component is expected is akeady 
an indication that the material in which they originate has been 
lost by the binary system, the time variability observed in coinci- 
dence with the eruption definitely testifies in favor of its circum- 
stellar nature. In this scenario, the blue-shifted features would be 
explained as arising in material swept by previous outbursts, as 
suggested for SN 2006X (Patat et al. 2007 ) following the mech- 
anism proposed by Wood-Vasey & Sokoloski ( 200 6 ), and hence 
buttressing the link between the progenitor system of that partic- 
ular SN and RS Oph. In this respect, although based on a small 
number of objects, the resemblance of RS Oph to a few Type la 
SNe is quite impressive (Fig. 10 1. 

Of course, many issues remain open, like the exact mech- 
anism that produces the weakening of the absorption features, 
the location and origin of the gas where the blue-shifted compo- 
nents arise (inhomogeneities in the RG wind vs. relics of previ- 
ous outbursts), and so on. But, as in the case of SNe, the study 
of strong absorption lines offers the unique possibility of investi- 
gating the presence of low amounts of neutral gas, which would 
never be detected through the typical signatures produced by di- 
rect interaction. Clearly, in the case of a Type la SN, for which 
ejecta velocities are much larger, the material generating com- 
ponent #3 would be swept during the first couple of days after 
the explosion, and hence it would be very hard (if not impossi- 
ble) to observe the same line variations we have detected in RS 
Oph (see the discussion in Williams et al. .2008. and Borkowski, 
Blondin & Reynolds 120081 about the possible detection of nar- 
row absorption features by means of very early, high-resolution 
spectroscopy of Type la. See also Hachisu, Kato & Nomoto l2008l 
for theoretical CSM scenarios in Type la progenitors). 

This is why well time-sampled, high-resolution spectroscopy 
of recurrent nova systems is crucial for better painting the com- 
plex picture that is just emerging. The recent report by Sternberg 
et al. ( 1201 11 1 of an excess of blue-shifted Nai absorption features 
along the lines-of-sight to Type la supernovae might indeed find 
a physical background within this context. 
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